Product quality evaluation by confidence intervals of process yield index

Statistical techniques have a beneficial effect on measuring process variability, analyzing the variability concerning product requirements, and eliminating the variability in product manufacturing. Process capability indices (PCIs) are not only easy to understand but also able to be directly employed by the manufacturing industry. The process yield index offers accurate measurement of the process yield, and it is a function of two unilateral six sigma quality indices. This paper initiates to develop the confidence intervals of the process yield index by using joint confidence regions of two unilateral six sigma quality indices for all quality characteristics of a product. Then integrate these joint confidence regions to find the confidence intervals of the product yield index. All manufacturing industries can use these confidence intervals to make statistical inferences to assess whether the process capability of the product and all quality characteristics has reached the required level, and to grasp the opportunities for improvement. An illustrated example on driver integrated circuit of micro hard disk is provided.

Process capability indices are commonly employed to assess whether the product quality can meet specifications defined in the manufacturing industry [1][2][3] . At the same time, the process capability indices are also commonly used in the industry, and many studies have invested in related discussions [4][5][6] . Based on some studies, a product usually contains multi-quality characteristics, including smaller-the-better (STB), larger-the-better (LTB), and nominal-the-best (NTB) at the same time 7,8 . Each quality characteristic needs to meet the required quality level, so that the quality of the final product can be guaranteed 9 . Moreover, numerous statisticians and quality engineers have studied process capability indices, aiming to come up with more effective methods to evaluate process potential and performance [10][11][12] . The six-sigma method is also a commonly used in the industry. Some studies are discussing the relationship between the six-sigma method quality level and the process capability index 13 , and then study and propose some six-sigma quality indices 14 . Two well-known unilateral six sigma quality indices, Q pu and Q pl proposed by Chang et al. 14 , are used to measure the STB and LTB quality characteristics as follows: and where USL is the upper specification limit, LSL is the lower specification limit, µ is the process mean, and σ refers to the process standard deviation. In normal condition, the process yield ( %Yield ) and unilateral six sigma quality indices Q pu and Q pl have close relations displayed below: and (1) where is a standard function of the normal cumulative distribution. The process yield and the unilateral six sigma quality indices own a one-to-one relationship in mathematics. Therefore, some studies use these two unilateral six sigma quality indices to develop fuzzy quality evaluation model 15,16 and fuzzy supplier selection model 17 .
For bilateral process capability index, Kane 5 proposed a yield built on capability index C pk below: Given that C pk is a function of µ and σ , it simultaneously depends on µ and σ . Hence, a confidence interval for C pk can be obtained using a joint confidence region for these two parameters. The lower confidence bound is the minimum value of C pk over the region. The approximate confidence bounds of C pk can then be obtained 18 . Hence, the process evaluation of C pk cannot provide exact process capability measurement and process yield.
Therefore, Boyles 19 has put forward a new bilateral process capability index which develops a one-to-one relation with the process yield as follows: Then, process yield index S pk offers accurate measurement of the process yield. When S pk = c , then the process yield is Yield = 2�(3c) − 1 . For processes of the normal distribution, the number of non-conformities is 2700 ppm, corresponding to a capable process with S pk = 1.0.
Index S pk has been widely used to judge whether the process quality can meet specifications and exactly measure the process yield. For example, Lee et al. 20 have proposed an asymptotic distribution for an estimator Ŝ pk . The asymptotic distribution of Ŝ pk is functional in statistical inferences for S pk . Huang et al. 8 have applied process yield index S pk to assess the product quality of a backlight module with multiple process characteristics. Whereas, first, the integrated product capability ( S T pk ) of the backlight module is defined and, second, the individual process quality for each quality characteristic is determined. In addition, Chen et al. 9 have considered generalizing the process yield index S pk targeted at processes with multiple quality characteristics. Wang and Du 21 used index S pk to assess the performance of the supply chain. According to index S pk , Wang et al. 22 developed a new index to assess the measurement of the yield for a multiple-stream process. Lin and Pearn 23 identified the problem of process selection with S pk so as to compare two processes as well as choose the one with a better production yield.
Unfortunately, the sampling distribution of yield index S pk is pretty complicate, and it is not easy to derive the confidence interval of the process yield index S pk . To overcome this difficulty, we have adopted several existing techniques to construct the confidence bounds for S pk [24][25][26][27][28][29][30][31] . For example, Chen 24 used the bootstrap simulation technique to figure out four approximate lower confidence limits of the yield index S pk . Shu and Wu 28 developed a useful method to gain the fuzzy estimate of the process yield index S pk for measuring the manufacturing process yield. Wu et al. 30 proposed a generalized confidence intervals for S pk to assess the process yield. However, most of the confidence bound studies for S pk are complicated and approximate estimation.
The yield index S pk is a function of indices Q pu and Q pl . Based on the Kushler and Hurley 18 method constructing the minimum of C pk over the region of the lower confidence bound, this paper initiates to develop the confidence regions of yield index S pk by using the joint confidence regions of the two indices Q pu and Q pl . This research aims to appraise the individual process quality of a multiple-process product according to the confidence intervals of the process yield index S pk . The method suggested by this study is a simple and intuitive tool. Then quality engineer can evaluate the process capability of the product and all quality characteristics, and decide whether to carry out process improvement. Therefore, the quality evaluation model in this paper can help the industry to improve.

Confidence intervals
Many studies have suggested that companies use control charts to perform process control. If the process is under statistical process control, then the process capability will be evaluated 32,33 . It is assumed that each subsample contains n observations on quality characteristics, and there are m subsamples available.
In each subsample, we let X i be the sample mean and S i be the sample variance of the i-th subsample, as displayed below: We define the overall sample mean and the pooled sample variance as follows: and The estimator Q pi of index Q pi is displayed as follows: For unilateral six sigma quality index Q pi , the 100 × (1 − α)% lower and upper confidence limits LQ pi and UQ pi for Q pi satisfy Based on Choi and Owen 34 , √ mn ×Q pi follows a non-central t distribution with m(n − 1) degrees of freedom, where n is the subsample size and m is the number of sample groups. The non-central parameter of Then, the above two equations can be rewritten as and where i = u or l and α ′ = α q ; q is the total number of quality characteristics for a multi-process product, m is the number of sample groups, and n is the sample size for each sub-sample. Online Appendix S1 displays the lower limit ( LQ pi ) and the upper limits ( UQ pi ) of the 95% confidence intervals for different Q pi values with n = 11 , m = 30 , and q = 6(1)10 . For instance, given sub-sample size n = 11 , m = 30 , q = 6 , α = 0.05 , when Q pi = 3.60 , from SAS program results, the LQ pi and UQ pi values are 3.222 and 3.981. The confidence intervals for different Q pi values can be obtained in Online Appendix S1. Online Appendix S1 only provides the numerical values which are used in this study. Different sample size, different number of sample groups or different number of quality characteristics could possibly occur in the practical application. For easier explanation, we only perform the selected sample size and selected number of sample groups and the results are listed in Online Appendix S1. The numbers of quality characteristics are computed from 6 to 10 in Online Appendix S1.
Assume there are q quality characteristics in a product, the process capability index for the jth characteristic will become Let the confidence intervals of indices Q puj and Q plj are denoted by LQ puj , UQ puj and LQ plj , UQ plj , respectively. And index S pkj is a function of indices Q puj and Q plj , the confidence intervals for index S pkj can be described as follows.
Lower confidence interval for the j th characteristic ( LS pkj ): Upper confidence interval for the j th characteristic ( US pkj ): Thus, the integrated process capability index for the entire product is Let the confidence intervals of indices S T pk are denoted by LS T pk , US T pk , and the confidence intervals for index S T pk can be described as follows. Lower confidence interval for the integrated product ( LS T pk ): As previously mentioned and based on Eqs. (16,17,19,20), these confidence intervals allow statistical inferences to be made to assess whether the product's process capability and all quality characteristics are at the required level.

Application procedures of confidence intervals
Products are commonly designed with many quality characteristics. Pearn et al. 25 has proposed how to determine the incorporated process of a product which contains multiple processes using index S pk . What we emphasize here is the construction and the application through the confidence intervals of the process yield index. The previous section illustrates how to calculate the confidence intervals for the incorporated process capability of a product which contains multiple processes on the basis of the process yield index S pk .
We assume that q quality characteristics in a product, m sample groups, and sub-sample sizes n are collected for each process. Since process yield index S pkj is a function of unilateral indices Q puj and Q plj , the estimates of Q puj and Q plj for each process should be calculated. Then the confidence intervals LQ puj , UQ puj and LQ plj , UQ plj for the estimates of Q puj and Q plj are determined according to the formulas in "Application procedures of confidence intervals". Or confidence intervals can also be found in Online Appendix S1. According to the confidence intervals for the estimates of Q puj and Q plj , the confidence intervals LS pkj , US pkj for each process is calculated. Finally, the confidence intervals for the incorporated process capability of the product with multiple processes are determined according to the formula of LS T pk and US T pk . For practical and easier application, the Fig. 1 evaluation procedure steps chart is shown as below.
The above evaluation steps are provided for manufacturers to follow and calculate the index confidence interval of the product and all quality characteristics, thereby evaluating whether the process capability of the product and all quality characteristics meets the quality level requirements. Next, this paper will use an example to illustrate the application of the above evaluation steps. Figure 1. Evaluation procedure steps chart.

An illustrative example
As the market demand on DVD player, MP3 and PDA rapidly grows, the demand on gigabyte capacity of multimedia data storage also increases. Nowadays, micro Hard Disk Driver (micro HDD) and Flash Memory card are the major products in terms of multimedia data storage. The storage capacity of a Flash Memory card ranges from 16 to 128 GB, while the storage capacity of micro HDD ranges from 500 GB to 5 TB. Driver IC is the key electronic component of micro HDD, functioning as anti-mechanical shock and read head control. In addition, it is responsible for the performance and effectiveness of micro HDD. The hard disk adopts the PCMCIA interface to connect to other equipment. To fit the design of thinner, smaller and low profile of micro HDD, the vertical dimension of driver IC gets thinner and thinner. There are six essential quality dimensions including top space, top loop, die, film, mold thickness and substrate thickness for a driver IC. The abovementioned evaluation method is adopted to assess the incorporated capability for a multi-process product, driver IC, for 1.8″ HDD with mold thickness of 0.25 mm. Table 1 shows brief information about this production line. We illustrate the above procedure step by step with a numerical example as follows: Step 1 Calculate the estimates Q puj and Q plj of indices Q puj and Q plj for each process in Table 2. For example, the Q puA and Q plA for process A are 2.73 and 4.08, respectively.
Step 2 Find the confidence intervals of indices Q puj and Q plj , which are denoted by LQ puj , UQ puj and LQ plj , UQ plj for each process in Table 2. For example, the LQ puA , UQ puA and LQ plA , UQ plA for process A are (2.430, 3.033) and (3.660, 4.503) , respectively.
Step 3 Calculate the confidence intervals LS pkj , US pkj for each process in Table 2. For example, the confidence intervals LS pkA , US pkA for process A are (2.667, 3.237).
Step 4 Calculate the confidence intervals for the integrated process capability for the 6-process driver IC. The confidence interval for the product is (2.888, 2.993).
Step 5 Compare the confidence intervals of driver IC with capable value, c = 1.0 . Since the capable value does not range in the confidence intervals, the integrated process capability for the driver IC is not capable.
Step 7 From Table 2, compare the required value C 0 with the confidence interval for each process and make a decision to determine which quality characteristic need to improvement. If C 0 ranges LS T pk and US T pk , then it is concluded that the individual process capability can meet the preset target; otherwise, the conclusion will be reverse. The processes marked with an " * " indicate that the process capabilities are not capable.
According to the above evaluation steps, the quality engineer can complete the process capability evaluation of the product and all quality characteristics, and decide whether to carry out process improvement. From Table 2, items B and F are capable processes. Quality engineers need to launch quality enhancement projects on incapable items A, C, D and E for process capability improvement. Obviously, through the evaluation steps in this  www.nature.com/scientificreports/ article, process engineers can simultaneously master and improve four quality characteristics with insufficient process capabilities. When the process capabilities of all quality characteristics meet the quality requirements, the product's process capabilities will meet the quality requirements.

Conclusions
PCIs are widely employed by the manufacturing industry to evaluate whether the process capability can meet the specifications. Regarding the product with multiple processes, customers concern the incorporated capability of the product. Based on the yield index S pkj for quality characteristic j, this paper discussed the incorporated process capability of a product with multiple processes in terms of calculating the confidence intervals of the yield index S pkj . Whereas the yield index S pkj is a function of indices Q puj and Q plj , the confidence intervals for individual process of indices Q puj and Q plj are computed to attain the confidence intervals of S pkj . Then integrate these confidence intervals to derive the confidence intervals of the entire product yield index S T pk . Evaluation procedures are presented in steps to assist practical application. The quality engineer can follow the evaluation procedures to complete the process capability evaluation of the product and all quality characteristics, and decide whether to carry out process improvement. The above research is based on the premise of the normal process. When the process distribution is non-normal, the method in this paper will have a large error, so it can be the focus of future research.